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The article is devoted to Be — >■ Bi*^ + mv decays with n = 1, 2, 3, 4. In the framework of 
factorization theorem the branching fractions of these processes can be written as convolution of 
hard part, describing Be Bi*''W vertices, and spectral functions, that correspond to transition 
of virtual W-boson into a final 7r-meson system. These functions were obtained from the fit of 
experimental data on r-lepton decay and electron-positron annihilation. Using different sets of 
Be Bi*'^W decay form-factors we present branching fractions and distributions over the invariant 
mass of the final 7r-meson system. 
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; I. INTRODUCTION 



i?c-nieson is the heaviest of the particles stable with respect to strong and electromagnetic interaction. The decays 
of the ground state of (6c) -system can be caused only by weak interaction, with either c-quark decays, 6-quark decays 
or annihilation allowed. According to [1], the ratios of these processes are ^ 45%, 37% and 18% respectively. Up to 
now only two decay modes, caused by 6-quark decay, are observed: Be — > J/ipir and semileptomic decay Be — J/ipiv. 
With the help of the former mode i?c-meson mass was determined with pretty good accuracy: niB^ = 6275.6±2.9±2.5 
^ , MeV (CDF collaboration [2|) and tob^ = 6300 ± 14 ± 5 MeV (DO collaboration [3]). The latter decay mode gives 
(]j ' the opportunity to measure i3c-meson life time: tb^ = 0.448lQg3g ± 0.032 ps (DO collaboration and tb^ — 
i 0.475^Q Q49 ^ 0.018 ps (CDF collaboration Q). Both the mass of _Bc-meson and its lifetime are in good agreement 

with theoretical predictions based on potential quark models and QCD sum rules [7l-[lT|. 
' According to perturbative QCD estimates |12l |. the cross section of Bc-nieson production in hadronic experiments 
J> ' is about 10"'^ of i3-meson production cross section. As a result, one can expect about 10^ Bc-mesons at LHC collider 
C ■ ', luminosity ~ lfb~^, so detailed investigation of ground i?c-meson and excited states of (6c)-family (there are 16 
narrow states below i?iD-pair production threshold) would be possible. 

In our previous work [13] we considered i?c-meson decays with 6-quark as a spectator, namely the reactions Be — > 
J/ip + X, where X stands for lepton pairs, light quarks ud or 7r-meson system m: for 1 < n < 4. In the framework of 
\l factorization theorem the amplitude of this process can be written as a convolution of Be J/ V'VF-decay width and 
spectral functions that describe the transition of virtual M^-boson into final 7r-mesons state. These spectral functions 
do not depend on H^-boson production mechanism, so they can be determined from experimental and theoretical 
. . , analysis of other reactions, for example r-lepton decay t ^ Vt + X ot electron-positron annihilation e+e" — ^ X. In 

. ^ ' the current work we consider induced by c-quark weak decay reactions Be — >■ i?s*'' + X. In contrast to Be J/^ + X 
: decays, these processes are Cabbibo-allowed, so their branching fractions are about an order higher than the branching 
\^ ' fractions oiBe^J/'(p + X decays. Form theoretical point of view these branching fractions are determined, from one 

5^ , side, by effective hamiltonian of weak c — s decay, where all higher-order QCD corrections are taken into account, 
and, on the other hand, by form-factors of Br — » Bl transitions. These from-factors can be obtained in different 
models, for example QCD sum rules H B-S [Hi , potential quark models [3, [S [13 , h ght-front approach fl9l-l23t. 
Calculated with these form factors widths Be ^ Bl + X decays differ from one another and, since considered 
here decays are dominant, this leads to different predictions for i?c-meson lifetime. As we have mentioned above, 
the lifetime of _Bc-meson is known experimentally with pretty good accuracy, so this difference can be used to select 
physical set of form-factors. 

In the next section we present analytical expressions for transferred momentum distributions of branching fractions 
of the decays Be ^ Be ' + X^ where X is lepton pair or light meson system. These distributions are expressed 
through spectral functions of the final state X. In section III explicit expressions of spectral functions, obtained 
from analysis of r-lepton decay r — > i^^ + X and electron-positron annihilation e+e^ — >■ are given. Using these 
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II. ANALYTICAL RESULTS 

Decays Be -B^*'' + mr are caused mainly by weak c-quark decay c — ?► sW~^ — > ud (see diagram shown in figHJ, 
while 6-quark is a spectator. Effective Lagrangian of this process has the form 



^eff = ^KsKd[C^+(M)0+ + C_(M)0-], 

where Gf is Fermi coupling constant, Vij are CKM mixing matrix elements and C±{fi) are Wilson coefficients, that 
describes higher-order QCD corrections. The operators 0±are defined according to 



0± = {d^UJ)^_^{c^bJ)y_^± {djU^)^_^{c,bj)y_^, 

where i, j are colour indexes of quarks and {qiq2)v-A ~ 9i7m(1 ~ 75)92- Since in the considered here decays the 
light-quark pair is in colour-singlet state, the amplitude of these processes should be proportional to factor 

«i(m) = ^[(A^c-l)C+(/i) + (7V, + l)C_(/i)]. 

If QCD corrections are neglected, this fuction is equal to ai (fj,) = 1. Due to higher-order logarithmic corrections the 
dependence of this coefficient on renormalisation scale fi appears j^l] , and on ^ ~ rUb it is equal to [To| 

0-1 (rriiy) = 1.2. 

The matrix element of the decay Be ^ Be +7?., where TZ is some set of light hadrons, has the form 



M 



B+ Bi*'>W+ b[*^ +n 



where is the effective polarization vector of light hadron system 7?,, and T-L vertex is 
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= .f+ (q^) Pf. + f- (q^) q^. 

for i?c-nieson decay to pseudo-scalar iJ^-meson and 



for Sc-meson decay to vector i3*-niesoii. In these expressions pi 2 are the momenta of Be- and Bs -mesons respectively, 
p = Pi +P2, q ^ Pi -~P2 is the momentum of virtual M^-boson, and f± {q^), Fif± (9^) and Fy (q^) are the form-factors 



of Be Bi*^W* transition. 



From presented above amplitudes it is easy to calculate the widths of the Be ^ Bs +TZ decays: 



dV 



Be Bi*'> + n 



1 GlV^, 



aln^n:e'j^e*T^d<i> (Be ^ B^*^ + U 



2M 2 

where M ~ Ms^is the mass of initial i?c-meson, and Lorentz-invariant phase space is defined according to 



(i$ (Q ^ fci . . . kn) = i2Trf S^I^Q-Y. j n 



(2^)' 2E, 



The following recurrence relation holds for this expression: 

d$ ^ B(*)7^) = ^d<^> (^Be ^ Bi*'^W*^ d^ {W* ^ TZ) , 
Using it one can perform the integration over the phase space of light hadrons system TZ: 

^jd^ {W* ^ n) e^e*J^ = {q^q, - q^g^,) (q^) + q^q^p^ [q^) . 

In the framework of factorization theorem introduced here spectral functions ^ (g^) are universal, so they can be 
determined from theoretical and experimental analysis of other reactions, for example r-lepton decay t — > Vt+TZ (25| 
or electron-positron annihilation e^e~ — >■ TZ. Explicit expressions of these spectral functions for different final states 
TZ are presented in the next section. 

From presented above matrix elements it is easy to obtain squared transferred momentum distributions for the 
considered in our article decays. In the case of pseudo-scalar Bs -meson we have 



dV{Be^ Bs+TZ) ^ GlYelaj 
32 



dq^ 



q^pf + 2Rc [/+/:] g2 (Af 2 - m') p^} , 



where M and m are the masses of Be- and Bs -mesons respectively, and 



' (M - mf - g2 / (A/ -I- mf - q^ 



M2 



M2 



is the velocity of B^-meson in i?c-meson rest frame. In the case of vector _Bs-meson in the final state the distribution 
has the form 
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SR [10] 


PM 


LF [26| 






SR [10] 


PM [10] 


LF [26] 




F(0) 


1.3 


1.1 


0.73 




F{0) 


8.1 


8.2 


6.1 


F+ 


ci, GeY-'^ 


0.30 


0.30 


0.56 


Fo^, GeV 


ci, GeV-2 


0.30 


0.52 


0.56 




C2, GeV"" 


0.069 


0.069 


0.030 




C2, GeV"'' 


0.069 


0.02 


0.087 




F(0) 


-5.8 


-5.9 


-1.7 




F{0) 


0.15 


0.30 


0.30 


F- 


ci, GeY-^ 


0.30 


0.30 


0.70 


F^, GeV"' 


ci, GeV-2 


0.30 


0.30 


0.30 




C2, GeV"* 


0.069 


0.069 


-0.02 




C2, GeV"'' 


0.069 


0.069 


0.069 




F(0), GeV"' 


1.08 


1.1 


0.31 




F(0) 


1.8 


1.4 


1.4 


Fv 


ci, 


0.30 


0.30 


0.12 


F^, GeV"' 


ci, GeV"2 


0.30 


0.30 


0.30 




C2, GeV"^ 


0.69 


0.069 


-0.02 




C2, GeV""* 


0.069 


0.069 


0.069 



Table I: Form-factor parameters 



dq^ " 1287r m^^^P^ 



2 „2\ o2t3„ /'C1AT;^A*^ 



+ Pl/?' [|Fo^|' + (Af2-m2)'|F^|' + g4|^A|2^ 

+ 2g2(A/2_^^2)j^g^^A^A*)]|_ 

One can easily see, the Fy (<z^) from- factor enter in this expression only quadratically, so it is not possible to determine 
its relative phase from transferred momentum distributions. 

There are several; theoretical models, that predict Be ~^ Bg + W decay form- factors: QCD sum rules [l^l, potential 
quark models [l^l, covariant light-front models [l^, etc. These from- factors can be parametrised in different forms, 
for example 

1. monopole expression 

^(0) 



Fdq' 



1-^2/^2' 



2. Isgur-Wise function £,{w) 



3. exponential parametrization, suitable for potential models, where form-factors are determined by integral of 
initial and final quarkonia wave-functions, that have Gaussian from: 

F,{q^) = F(0)exp{ci<72 + C2g4}. (1) 

All parametrizations are almost equivalent and in our article we use this the exponential parametrization ([T]) for all 
sets of form-factors. Numerical values of the parameters F (0), ci,2 for different from-factors stets are given in tableU 

III. NUMERICAL RESULTS 

When inclusive decays Be Bg ud are consideres, corresponding spectral fiuctions are equal to 

pr il') = 0, p-r' id') = (2) 
Branching fracttions of these decays for listed above form-factor sets are 



BrsR {B,ud) = 19%, BrpM (Bsud) = 13%, BrLF (Bsud) = 6.5%, 
BrsR {B*ud) = 20%, BrpM {B*ud) = 23%, BrLF {B*ud) = 13%. 
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The braching fractions of semileptonic decays Be — >■ -Bs eVe can be obtained from these values by simple substitution 

Let us now consider Be Bs decay. The IF — )■ tt vertex is written as 

(vr+|J^|W^) = /.q^, (3) 

where is vr-meson momentum and constant can be determined, for example, from the width of leptonic decay 
TT"*" — e^v + e. Using experimental value for the width of this decay we obtain /tt ~ 130 MeV. Spectral functions, 
that correspond to vertex ([3]) are 

Pi {q') = fj {q' - ml) , (q^) = 0. 

The values oi Be ^ Bs n branching fractions for listed in the previous section form- factor sets are 

BrsR,(B,7r) = 18%, BrpM (B.tt) = 12%, BrLF (S.tt) = 5.5%, 
BrsR(S» = 7%, BrpM (B^tt) = 9.4%, BrLF (S» = 6.2%. 

It is clearly seen, that the difference in form-factors leads difference in the branching fractions of these decays. All 
these branching fractions, in the other hand, are about an order of magnitude higher, that the branching fractions of 
Be J/V'TT decay [l^. The reason is that on the quark level Be — J> -Bs*' + Tl decays are caused by Cabbibo-allowed 
c — sud decay, while Be J/ip + 71 decays are Cabbibo-suppressed. 

If there are two 7r-mesons in the final state, the main decay mode would be Be B^*^ p'^ B^f^^n^ . The vertex 
of p-meson interaction with VF-boson has the form 

where rrip and e^are p-meson mass and polarization vector, while /p-constant is fp « 210 MeV. In the limit of zero 
p-meson widths the spectral functions p^^^ (q^) are equal to 

p£ (q^) = 0, p?, (q^) = fl5 {q^ m^) . 
These spectral functions lead to following branching fractions oi Be ^ Bl p decays: 



BrsR [Bsp) = 7.6%, BrpM (-B^p) - 5.4%, BrLF [BsP) = 3.1%, (4) 
BrsR {BIp) = 21%, BrpM (B^) = 22%, BrLF {Blp) = 15%. (5) 

As in previous case, these values exceed significantly the branching fractions of Be — !■ J /i^p decay. 

It should be noted, however, that in contrast to i?c J/i^ + In decay, for Be ^ Bs + 2Tr decays it is not valid 
to neglect the width of p-meson. Because of the small mass difference Mb^ — Mb^ ~ 1 GeV p-meson is almost at the 
end of the phase space, so its non-zero width changes significantly g^-distributions for these decays. For this reason 
one should use more realistic parametrization for spectral functions pj^x ^'^^ to vector current conservation the 
longitudinal spectral function p|^ (g^) = 0. The information on transverse spectral function can be obtained from 
analysis of r-lepton decay t ^ + 27r. For this decay q^-distribution has the form 



dV (r ^ ^^^7^) G% (m^ - q 



2\ 



dq^ IGnmr 



F \ r ( 2 



2q') P? {q') , (6) 



where rUr is r-lepton mass, and, in the framework of factorization theorem, spectral function p|r (q^) is universal, 
that is independent on 7r-pair production dynamics. The authors of paper [27| give the parametrization for this 
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spectral function with p-, p'- and cj-meson contributions taken into account. In our article we use a more simple 
parametrization 



p^^ (s) = 1.35 X 10 



-3 



4m^ 



1 + 0.64s 



(s- 0.57) +0.013 



where s is measured in GeV^ (the spectral function itself is dimensionless) 

Distributions of — )■ si*^ + 27r-decays branching fractions over squared transferred momentum for listed in the 
previous section form-factor sets are presented in figH] Solid, dashed and dash-dotted lines in this figure correspond 
to form-factor sets "SR", "QM" and "LF" respectively. Corresponding branching fractions are 



BrsR (B.-KTi) = 6.1%, BrpM (S^tttt) = 4.3%, BrLF (BsTttt) = 2.4%, 
BrsR (B^tttt) = 13%, BrpM (S^tttt) = 14%, BrLF (S^vrvr) = 8.3%. 



(7) 

(8) 



One can easily see that these branching fractions are smaller than presented above values ([S]). In fig|3]we present 

q^-distributions of B^. — )■ B'f^ +271 decay (solid line) and Be ~^ iJs*"* +ud (dashed line. In this case the spectral function 
does not depend on q^, see eq.Q ). Experimental value of p- meson mass is shown on this figure by vertical lines. 
From this figure it is clear, that in the region « (mp zLT , where spectral function is significantly non-zero, 
the distributions dT (^Bc — > B^*'' + udj /dq^ vary strongly, so one cannot neglect p-meson width in this case. This is 

the reason for large difference between ([4]), ([5]) and (O, ([8]) branching fractions. 

If there are three 7r-mesons in the final state the main production mode would be ai — > p7r — > 37r. There are 
two possible charge configurations (7r+7r^7r+ and 7r+7r°7r°), and in our article 37r stands for sum of these final states. 
Because of p artial conservation of axial current longitudinal spectral function p^{q^) could be set equal to zero. 
In paper [27| the parametrization of transverse spectral function p"^ (g'^) , expressed through mass and width of oi- 
meson is proposed. In our article we use a simpler parametrisation, obtained from fit of experimental q^ distribution 
of T — >■ i/T- -|- 37r decay: 



p^"" (s) = 5.86 X 10" 



where s is measured in GeV . Distributions of Br 



9mi 



1 + 190s 



B. 



(*) 



(s- 1.06)V0.48 



37r branching fractions over q^ for different stets of 



form-factors are shown in fig|31 Corresponding branching fractions are equal to 

BrsR {B,3tt) = 0.096%, BrpM (BsSn) = 0.068%, BrLF (S,37r) = 0.039%, 
BrsR iB*37r) ^ 0.23%, BrpM (B^Stt) = 0.24%, BrLF iB*^Tr) = 0.16%. 

One can clearly seen, that these values are significantly smaller, than presented above branching fractions of the 



decays Be — M*"* -I- tt, B, 



B 



(*) 



2tt. This decays, however, could be interesting from experimental point of view. 
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Figure 3: Distributions of branching fractions of Be ^ Ba + 2tt, Be Bs + ud decays (left panel) and Be ^ Bg + 2tt, 
Be ^ B* + ud (right panel) over the squared transferred momentum. On both figures solid and dashed lines line correspond 
to TTTT and ud final states respectively. Experimental value of p-meson mass is shown by vertical lines. 




Figure 4: Distributions of Be — >■ Bs* + Stt branching fractions over squared momentum of 37r-meson system. Notations are 
same as in fig [2] 



since in tt+tt 7r+ charge configuration 7r°-meson, whose registration could be problematic, is absent. Moreover, in 
contrast to Be — > B^*^ + Be ^ B^*^ + 27r decays, the branching fractions of Be — >■ si*'' + Stt decay are smaller 
than the corresponding branching fractions of Be — >• J/ip + Stt decay. The reason is that for Bs -meson in the final 
sate the masses of initial and final mesons are rather close, so phase-space suppression compensates the enhancement 
caused by CKM matrix element. 

If there are four 7r-mesons in the final state two charge configurations are possible: tt^tt^tt^tt'^ and tt^ tt'^ tt^ tt^ , in 
what follows Att stands for sum of these configurations. Due to vector colour conservation the longitudinal spectral 
function Pj^ {q^) = 0. Information about traverse spectral function (^q^) is available from experimental data on 
T ^ + 4:TT decay (see eq.([n]) ) or electron-positron annihilation e+e^ — > Att. Cross section of the latter reaction is 

s 

In both cases we have similar from of the spectral function, that can be parametrized by the expression 

s - 16ml \ 1 - 5.07s + 8.635^ 
^ ^ [(s - 1.83)^ + 0.6lj^ 

In figlUg^-distributions of Be — > B^*^ +4n decay branching fractions for different form-factor sets are shown. Numerical 



1.8 X 10" 



dBr 1 rfBr 1 

dQ- GeV^ dQ^ GeV- 




Figure 5: Distributions of Be — >■ -Bs + 47r branching fractions over squared momentum of 47r-meson system. Notations are 
same as in figO U 



Br(Bc ^Be + 7^), % 


Br (Be ^B*+7^), % 


7^ 


SR [10] 


PM [10] 


LF [26] 


7^ 


SR [10] 


PM [10] 


LF [26] 


TV 


18 


12 


5.5 




7 


9.4 


6.2 


P 


7.6 


5.4 


3.1 


P 


21 


22 


15 


2ti 


6.1 


4.3 


2.4 


2tt 


13 


14 


8.3 


Stt 


0.096 


0.068 


0.039 


3ir 


0.23 


0.24 


0.16 


47r 


0.0064 


0.0045 


0.0026 


A-K 


0.015 


0.016 


0.01 


ud 


19 


13 


6.5 


ud 


20 


23 


13 



Table II: Branching fractions of Be — >■ -Bs + 7?. decays 



values of these branching fractions are 

BrsR (B,47r) = 0.0064%, BrpM (S,47r) = 0.0045%, BrLF [BAt^) = 0.0026%, 
BrsR (S*47r) = 0.015%, BrpM {B*At^) ^ 0.016%, Bilf (BI47r) = 0.01%. 

These values are significantly smaller then presented above branching fractions of one, two, or three Tr-meson pro- 
duction, as well as branching fractions of i?c >//V' + 47r decay [13]. The reason is mentioned above phase-space 
suppression m Be B'-*'^ + in reaction. Thus we can see, that for all form-factor sets dominant decay mode c — > s 
is saturated by Be ^ Bl n, p decays. 

In table |TT] branching obtained in our article fractions of B^ — ?> b[*^ -|- TZ decays for different form-factor sets are 
listed. It is interesting to compare these results with simple estimates based on duality relations. According to paper 
[H Bc-meson decays with spectator 6-quark takes about 45% of total _Bc-meson width. These decays are almost 

saturated hy B^ ^ Bs + X decays, so the sum oi Be ~> Bl + branching fractions should be about this value. 
It can be easily seen, that this is true for "SR" and "PM" form-factor sets (the corresponding sum for them is 44% 
and 40% respectively). For "LF" form-factors set this sum is significantly smaller 23%). From table [ll] it can be 
seen also, that for every decay mode the branching fraction obtained with "LF" form-factors sets is smaller then those 
obtained with "SM" and "PM" form-factors sets. As a result, the total width of B^ — > B^*^ + X inclusive decay (and, 
hence, lifetime of _Bc-meson) would contradict experimental value. 



IV. CONCLUSION 



The article is devoted to study of exclusive Sc-meson decays with production of S^* -meson an light hadron system 
rm with n = 1, 2, 3 or 4. 

i?c-mesons, that is particles that in valence approximation are build from c- and &-quarks, take the intermediate 
place between charmonia (cc-mesons) and bottomonia (66-mesons), so they can be used for independent check of 
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theoretical models used for analysis of heavy quarkonia with hidden flavour. Only ground state of i?c-nieson family 
is now obsevrved experimentally, with only mass and lifetime measured with good accuracy. One can expect large 
_Bc-meson production at LHC collider, so branching fractions of exclusive i?c-meson decays could be measured and 
theoretical prediction are desirable. 

Previous works were mainly devoted to two-particle decays of _Bc-meson (see, for example, 0, S H, [ll, [H-ilj ) . In 
our recent article [isj we, on the contrary, consider decays Be — > J/ip + with n = 1, 2, 3 or 4. In the framework 
of factorization theorem the branching fractions of these decays are written as convolution of hard part, describing 
Be — > J/ipW decay and spectral functions, responsible for W — nir transition. Form factors of Be — J/tpW vertex 
can be determined from different theoretical models (QCD sum rules, potential models, light-front covariant quark 
models, etc), spectral functions — from analysis of r-lepton decays and electron-positron annihilation. 

In the present article we consider in the same approach the decays B^. ^ Bl ' +mT with n = 1, . . . , 4. Using different 
sets of Be — > Bi*^W vertex form-factors we calculated branching fractions of these decays and distributions over the 
squared momentum of 7r-mesons system. In contrast to Be — >■ J/tp + nn decays these reactions are Cabbibo-allowed, 
so for one and two 7r-mesons in the final state their branching fractions are greater, than the branching fractions of 
corresponding Be — > J/ip + mr decays. If the number of 7r-mesons in the final state is larger, the suppression caused 

by small phase-space in the reaction Be Bs + X is important and Be ^ J/ip + X decays dominate. 

We would like to note, that considered in our article decays are suitable for i3c-meson study at hadron colliders. 
Besides large branching fractions there are also small background processes to these decays. 

This work was supported in part by Russian foundation for basic research (grant #10-02-00061a). One of the authors 
(A.V.L.) was also supported by Russian Foundation of Science Support, non-commercial foundations "Dynasty" and 
the grant of the president of Russian Federation (grant #MK-406.2010.2). 
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